Chemical zonation in olivine-hosted melt inclusions

Significant zonation in major, minor, trace, and volatile elements has been documented in naturally glassy olivine-hosted melt inclusions from the Siqueiros Fracture Zone and the Galapagos Islands. Components with a higher concentration in the host olivine than in the melt (e.g., MgO, FeO, Cr2O3, and MnO) are depleted at the edges of the zoned melt inclusions relative to their centers, whereas except for CaO, H2O, and F, components with a lower concentration in the host olivine than in the melt (e.g., Al2O3, SiO2, Na2O, K2O, TiO2, S, and Cl) are enriched near the melt inclusion edges. This zonation is due to formation of an olivine-depleted boundary layer in the adjacent melt in response to cooling and crystallization of olivine on the walls of the melt inclusions, concurrent with diffusive propagation of the boundary layer toward the inclusion center. Concentration profiles of some components in the melt inclusions exhibit multicomponent diffusion effects such as uphill diffusion (CaO, FeO) or slowing of the diffusion of typically rapidly diffusing components (Na2O, K2O) by coupling to slow diffusing components such as SiO2 and Al2O3. Concentrations of H2O and F decrease toward the edges of some of the Siqueiros melt inclusions, suggesting either that these components have been lost from the inclusions into the host olivine late in their cooling histories and/or that these components are exhibiting multicomponent diffusion effects. A model has been developed of the time-dependent evolution of MgO concentration profiles in melt inclusions due to simultaneous depletion of MgO at the inclusion walls due to olivine growth and diffusion of MgO in the melt inclusions in response to this depletion. Observed concentration profiles were fit to this model to constrain their thermal histories. Cooling rates determined by a single-stage linear cooling model are 150–13,000 C h from the liquidus down to *1,000 C, consistent with previously determined cooling rates for basaltic glasses; compositional trends with melt inclusion size observed in the Siqueiros melt inclusions are described well by this simple single-stage linear cooling model. Despite the overall success of the modeling of MgO concentration profiles using a single-stage cooling history, MgO concentration profiles in some melt inclusions are better fit by a two-stage cooling history with a slower-cooling first stage followed by a faster-cooling second stage; the inferred total duration of cooling from the liquidus down to *1,000 C ranges from 40 s to just over 1 h. Based on our Communicated by Jon Blundy. Electronic supplementary material The online version of this article (doi:10.1007/s00410-014-1030-6) contains supplementary material, which is available to authorized users. M. E. Newcombe (&) A. Fabbrizio C. Ma M. Le Voyer Y. Guan J. M. Eiler E. M. Stolper Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, CA 91125, USA e-mail: [email protected] A. Fabbrizio Laboratoire Magmas et Volcans, CNRS UMR 6524, Université Blaise-Pascal, OPGC-IRD, 5 Rue Kessler, 63038 Clermont-Ferrand Cedex, France Y. Zhang Department of Earth and Environmental Sciences, University of Michigan, Ann Arbor, MI 48109, USA M. Le Voyer Department of Terrestrial Magnetism, Carnegie Institution of Washington, Washington, DC 20015, USA A. E. Saal Department of Geological Sciences, Brown University, Providence, RI 02912, USA 123 Contrib Mineral Petrol (2014) 168:1030 DOI 10.1007/s00410-014-1030-6